Contrast of LiFeAs with isostructural, isoelectronic, 
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Stoichiometric LiFeAs at ambient pressure is an 18 K superconductor while isoelectronic MgFeGe 
is not, despite their extremely similar electronic structures. To investigate possible sources of this 
distinctively different superconducting behavior, we quantify the differences using first principles 
density functional theory, establishing first that the Fe total 3d occupations are identical in the two 
compounds. Individual 3d orbital occupations also differ very little (~ 0.01). The differences in 
Fermi surfaces (FSs) do not seem significant; however a redistribution of bands just above the Fermi 
level does represent a possibly significant distinction. Because the bands and FSs of LiFeAs are 
less in agreement with experiment than for other iron-pnictides, we study the effects of additional 
exchange-correlations effects beyond GGA (the generalized gradient approximation) by applying the 
modified Becke- Johnson potential (mBJ) exchange potential, which gives much improved bandgaps 
in insulators compared to GGA and might be useful for semimetals such as the Fe-based supercon- 
ductors. Overall, we conclude that the mBJ corrections do not improve the description of LiFeAs 
as compared to experiment. 

PACS numbers: 71.20.-b, 74.20.Pq, 74.70.Xa 



I. INTRODUCTION 

In the four years since the discovery of super- 
conductivity in the iron-pnictidc and -chalcogenide 
superconductors (FeSCs), a great deal has been 
learned about the materials physics of the handful of 
structural subclasses that comprise these new high 
temperature superconductors. Differences between 
the subclasses have been uncovered, but agreement 
on the microscopic mechanism of pairing is lacking. 
We suppose (as is commonly held) that the super- 
conductivity that occurs in these materials which 
have in common a layered, fluorite-type Fe-Pn back- 
bone (Pn = pnictide or chalcogenide) has a com- 
mon origin, at least for the "high T c " (> 10-15 K) 
cases. One of the greatest current needs is to iden- 
tify microscopic characteristics that can shed light 
specifically on the existence, or not, of high temper- 
ature superconductivity (HTS) and thereby on the 
underlying pairing mechanism. 

Several mechanisms have been put forward. Be- 
cause superconductivity borders and competes with 
magnetic order as in the HTS cuprates, it is natural 
to study a spin fluctuation (SF) origin, and several 
groups^— have pursued antiferromagnetic SF mod- 
els. A ferromagnetic SF model has been suggested 
by Brydon et al.^ largely on the basis on the poor 
nesting of its calculated FSs. Fe 3d orbital occu- 
pation and character have received much attention, 
and an orbital fluctuation (OF) model has been sug- 
gested by Saito et al& The role of the Pn anion 
was given more attention in the charge fluctuation 
(CF) picture of Zhou et al.^ where interatomic Fc- 
Pn charge-charge interactions provided another elec- 



tronic mechanism. On one hand, the mixing of the 
narrow Fe 3d bands removes almost all Pn charac- 
ter in the states near the Fermi level sp, making 
most models focus simply on the 3d states. On the 
other hand, T c has been found to correlate strongly 
with the distance of the Pn above and below the 
Fe plane, or more specifically on the Pn-Fc-Pn bond 
angle and/or Fe-Pn bond length. 

The Pn anion received attention early in the 
study of FeSCs, when Berciu et cd4 modeled pair- 
ing in terms of electronic polarons and bipolarons. 
Their picture foreshadowed the CF model mentioned 
above, but placed more emphasis on the anticipated 
large polarizability of the Pn anions, and did not 
include consideration of SFs or OFs. This picture 
was continued and extended to include interatomic 
exchange coupling by Chan et al£ The conventional 
phonon mechanism, for which there is a reliable mi- 
croscopic theory when electronic interactions are de- 
scribed sufficiently by density functional methods, 
has been evaluated by Boeri et al^ and found to be 
too weak to explain T c in the 25-55 K range. Strong 
spin-lattice coupling, however, has been suggested 
by Egami et alJ^ to be involved in pairing. The 
electronic structure of these FeSCs is being studied 
in detail experimentally and modeled carefully by 
numerous groups, with a primary aim being to iden- 
tify the pairing mechanism. 

LiFeAs is recognized as a problem child in the cat- 
egorization of FeSCs. Among the vast collection of 
FeSCs that have been discovered since the first of 
their kind^ LiFeAs is one of very few which super- 
conduct without the need for either chemical dop- 
ing or physical compression. Of these few, LiFeAs 



not only has the highest superconducting transition 
temperature (T c = 18 K), but it is to date the only 
compound, other than LiFeP, whose T c does not in- 
crease when doped or pressurized. It is widely be- 
lieved that in the majority of FcSCs, doping or ap- 
plication of pressure suppresses the nesting-induced 
spin-density-wave order to make way for compet- 
ing superconductivity. LiFeAs however does not un- 
dergo any magnetic transition, and significance of FS 
nesting in LiFeAs has been questioned^ Not only 
does LiFeAs not fit well with many aspects of FcSCs, 
its behavior also differs strongly from its isova- 
lent sister compound NaFcAs. Similarities between 
NaFeAs and LiFeAs in band structure and den- 
sity of states are even more profoun d 14 ^ 15 than be- 
tween MgFeGe and LiFeAs, but unlike MgFeGe and 
LiFeAs, NaFeAs undergoes a magnetically driven 
structural phase transition above the superconduct- 
ing transition J£~— In fact, bulk superconductivity 
may not even exist in NaFeAs J£r— 

The recent report of the synthesis and char- 
acterization by Hosono's group^ 2 - of non-magnetic 
MgFeGe, which is isostructural and isoelectronic 
with the 18 K superconductor LiFeAs but is not 
superconducting, provides a means to obtain new 
insight. In their initial report, Liu et al£^ noted 
the resemblance of the electronic structure near ep 
to that of LiFeAs, providing both a conundrum and 
an opportunity to identify differences that account 
for the vast distinction in superconducting behav- 
ior. Supposing that Mg gives up both of its valence 
electrons to the Fe-Ge bands, two initial basic ques- 
tions emerge: the 3d charge on Fe, and the distinc- 
tions between the anions As and Ge that are neigh- 
bors in the periodic table. We first establish that 
the Fe 3d occupation is identical in these two com- 
pounds. The (presumably) more highly charged Ge 
anion should have an even higher polarizability than 
As, for which the simplest viewpoint might suggest 
to be more favorable for superconductivity, rather 
than precluding it. Moreover, the As-Fe-As angle, 
which has been shown for most of the pnictide o 23 ' 24 
to correlate strongly with increased T c as it ap- 
proaches the regular tetrahedral angle of 109.47°, is 
103.1°, very similar to MgFeGe's Ge- Fe-Ge angle of 
103.6°. (Note that NaFeAs, which possesses an angle 
of 108.3°, is a departure from these two compounds 
in this respect.) One can also wonder whether the 
alkali/ alkali earth atom can make any real difference. 
Otherwise, the essential difference would come down 
to small differences on the Fe atom, such as individ- 
ual orbital occupations, which, as we show, are quite 
small. 

In this paper we perform a close comparison of 
the electronic structures of MgFeGe and LiFeAs. 
The differences of the Fe 3d orbital occupations are 



small but clear. Because the local density approx- 
imation (LDA) and generalized gradient approxi- 
mation (GGA) have deficiencies in their descrip- 
tion of the magnetism and atomic structure of these 
FeSCs, generally and LiFeAs in particular, we as- 
sess the effect of using the modified Becke-Johnson 
(mBJ) form of exchange-correlation potential 2 ^ that 
builds in "exact exchange" in a local-density man- 
ner. The mBJ potential has now been established to 
provide excellent bandgaps in many semiconductors 
and insulators. Although it is not clear from most 
band structures, the FeSCs are close to a semimctal- 
lic state. Lcbcguc et al£&, in their Fig. 15, show 
that in the as yet unsynthesized LaOFeN the occu- 
pied and unoccupied bands form disjoint complexes. 
This observation raises the question whether these 
same exact-exchange corrections might apply to a 
semimctal. 



II. ELECTRONIC STRUCTURE 
A. Description of methods 

The non-magnetic electronic structures of LiFeAs 
and MgFeGe were calculated using WIEN2K, 2 ii 2 ^ 
a density functional theory (DFT)-based full poten- 
tial, linearized augmented planewave (LAPW) code. 
For the exchange-correlation functional, the gener- 
alized gradient approximation (GGA) of Perdew, 
Burke, and Ernzerho^ was applied. A fc-point mesh 
of 32 x 32 x 19 was used, and RK max was set to 10. 
As for the muffin-tin (MT) radii R, we used val- 
ues of 2.50, 2.43, and 2.14 a.u. for Li/Mg, Fe, and 
As/Ge, respectively. We adopted the experimental 
lattice constants of LiFeAs and MgFeGe (both of 
which form a PA/nmm tetragonal unit cell) given in 
Refs. and HO respectively. The internal coordi- 
nates came from the same experiments, and relax- 
ation was not applied in our calculations since we are 
interested in identifying distinctions and comparing 
with experiment when possible. The structural pa- 
rameters are given in Table |U Note that the Li (Mg) 
atom is five-coordinate with As (respectively, Ge) in 
a square-based pyramid, and its height above the Fe 
plane is 2.20 A (2.13 A). 

B. Bands and densities of states 

The densities of states (DOSs) of LiFeAs and 
MgFeGe are displayed in Fig. [TJ In the case of 
LiFeAs, the Fermi energy ep sits on the shoulder of 
a high DOS region with almost entirely Fe d char- 
acter, consistent with results presented by Singh,— 
Nekrasov et al.^ and Liu et al~ Li/Mg contribu- 



2 




Energy (eV) 



FIG. 1: (Color online) Densities of states (per two spins, 
per primitive cell) of LiFeAs and MgFeGe. Li/Mg states 
give negligible contribution in this energy range so they 
have not been plotted. The primary difference occurs in 
the 0.0-0.5 eV regions (see text). 

tions are negligible at and around the Fermi energy, 
and As / Ge density is also very low at ef ■ The large 
DOS at ef, N(ef), is a distinctive feature of LiFeAs 
compared to the smaller values in the other iron- 
pnictides, and MgFeGe shares this large value. In 
fact, N(ef) is 25% larger in MgFeGe than in LiFeAs, 
a property that should enhance superconductivity 
independent of mechanism rather than hinder it. 
While a set of experimental data has been used to 
characterize LiFeAs as intrinsically overdoped^ to 
explain its 18 K critical temperature, from its larger 
N(ef) value due to Ef falling further below the val- 



TABLE I: Crystal structure parameters of LiFeAs and 
MgFeGe in the tetragonal space group P4/nmm. Wyck- 
off labels of Li/Mg, Fe, and As/Ge are, respectively, 2c, 
2b, and 2c. 



a (A) 


C (A) 








Li 


(1/4, 1/4,0.8459) 


LiFeAs 3.7914 


6.3639 


Fe 


(3/4, 1/4, 1/2) 






As 


(1/4, 1/4,0.2635) 






Mg 


(1/4, 1/4,0.8316) 


MgFeGe 3.8848 


6.4247 


Fe 


(3/4, 1/4, 1/2) 






Ge 


(1/4, 1/4,0.2620) 



ley in the DOS it would seem to be underdoped. 

The LiFeAs and MgFeGe DOSs share similar Fe d 
features from —2 eV up to Ep , at which point, in 
MgFeGe, the shoulder is more abruptly cut off. The 
subband DOS peak of unoccupied states centered at 
0.5 eV in LiFeAs is shifted down in MgFeGe, to form 
a denser manifold just above the abridged shoul- 
der, and is centered 0.2 eV above the Fermi level. 
The separation of the Fe DOS into orbital compo- 
nents will be given later, when effects of different 
exchange-correlation functional are discussed. 

The band structures of the compounds are com- 
pared in Fig. 21 As noted by Liu et al., the sim- 
ilarity of their bands near ef is the most striking 
feature. The doubly degenerate band at T drops by 
less than 0.1 eV to lie exactly at the Fermi level in 
MgFeGe, and the pocket at M shrinks in size due 
to a rise by less than 0.1 eV of the band just be- 
low ef] the bands elsewhere at Ef are extremely 
similar. The increased concentration of states in 
MgFeGe just above the Fermi level, visible in the 
DOS as mentioned above, can be traced to the low- 
ering by 0.38 eV (from 0.43 eV to 0.05 eV) of the 
lowest unoccupied band at the X point. The FSs 
of LiFeAs and MgFeGe, discussed in more detail in 
later sections and shown in Fig. [7] and Fig. [Sj are 
each made up of three hole pockets centered at T 
and two electron pockets at M. The electron FSs and 
the largest hole pocket increase in size in MgFeGe, 
while the two smaller electron pockets shrink. We 
have verified that changes due to spin-orbit coupling 
are uninterestingly tiny. 



C. Charge density differences 

The Fe 3d occupations in LiFeAs and MgFeGe 
are identical. This essential feature is established by 
comparing the radial valence charge density r 2 p(r) 
of Fe between the two compounds, shown in Fig. [3] 
In the region of the 3d peak at radius r = 0.7 a.u., 
where the density is due only to 3c? occupation and 
to identical core tails, the difference is incredibly 
small (less than 0.1%), and the difference remains 
tiny out to 1.4 a.u. Using the same MT radius of 
2.43 a.u., the charge of Fe in this neighborhood in 
LiFeAs is larger only by 0.01 electrons, with this 
tiny difference arising from the differing As and Ge 
tail charges. The orbital occupation matrix elements 
fimm' for both compounds are listed for comparison 
in Table ILT1 No difference is larger than 0.015, which 
is the change for the d z 2 (m = 0) orbital, being 
smaller in MgFeGe. This difference is compensated 
by the d xz / yz occupations (m = ±1) increase by 
~ 0.01 in MgFeGe. The m = ±2 (d xy , d x %- y %) or- 
bital occupation changes are negligible. The t^g and 
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e g degeneracies are broken by the Fe site symmetry, 
of course, as can be noted in the occupation matri- 
ces. 



D. As- and Ge-related differences 

Since the total Fe 3c? occupation and even individ- 
ual orbital occupations differ little between the two 
compounds, it becomes of more interest to compare 
differences ascribable to As and Ge. One difference 
is the structure itself: the lattice constants a and c of 
MgFeGe are about 2.5% and 1% larger respectively. 
Of Fe's three nearest- neighbor interatomic distances, 
the Fe-Fe distance changes most, from 2.68 A to 
2.75 A, the 2.5% increase directly related to the 
same relative increase in a lattice constant. Effec- 
tive d-d hopping amplitudes might be changed some- 
what, however, since the hopping is largely through 
the As or Ge atom and this interaction will differ. 
There is a similar 2% increase in the Fe-As/Ge dis- 
tance (2.42 A/2.47 A), related to Gc's slightly larger 
atomic radius. Despite MgFcGc's larger unit cell 
(6% in volume), and Mg having larger nuclear charge 
than Li, but perhaps because Mg has an additional 
valence electron, the Fe-Mg distance is 0.8% shorter 
than that of Fe-Li (2.91 A versus 2.88 A). 

The Pn-Fe-Pn angles have been observed to cor- 
relate very strongly with the maximum value of T c 
observed within a structural class while varying sto- 
ichiometry and also across structural classes ! 23 ' 24 
This correlation does not extend well to the Fe 
chalcogenides, however, so its relevance in this com- 
parison of LiFeAs to MgFeGe may be questioned. 
The similarity of the angles for these two compounds 
seems to make the bond angle a non-issue in the 
question of presence or absence of superconductiv- 
ity. 

There has been interest in the polarizability of the 
pnictide atom and its effect on carriers. The polariz- 
ability of the Ge atom is 40% larger than that of As 
(6.07 vs. 4.31, in 10 -24 cm 3 ). Naively, this difference 
would favor superconductivity in MgFeGe according 
to the relevant models Moreover, with the addi- 
tional valence electron (Mg vs. Li) not going to Fe 
and thus going primarily to Ge, this difference in po- 
larizabilities of Ge and As may be enhanced in the 
solid, making it even more difficult to understand 
the difference in superconducting behavior. 



E. Susceptibility Xo(c\) 

In Fig. [2J we show the calculated constant- matrix- 
element bare susceptibilities Xo(*l) (normalized by 
Xo(q = 0)) of LiFeAs and MgFeGe, accounting for 



a) LiFeAs 




FIG. 2: (Color online) Non-interacting spin susceptibil- 
ities Xo{lx,q y ,qz = 0), each normalized by the xo value 
at q = (0,0,0), of (a) LiFeAs and (b)MgFeGe in the 
complete Brillouin zone. In these plots, the F point is 
positioned at the corners. 



all bands that cross the Fermi level. Our suscep- 
tibility for LiFeAs is consistent with that of Lee et 
al.j& who also observe a maximum in \o at (ir,ir). 
The two F-centered cylindrical FSs with radii fc_F,i, 
(see figures below) lead to the circular hump in 
the range of 2fcp i x, 2fcp,2, and /cf.i + Small 
differences in the band structure broaden the (ir, tt) 
peak in MgFeGe, but the overall xo topology is sim- 
ilar to that of LiFeAs. Applying dynamical mean 
field theory (DMFT) to DFT, in which dynamical 
electron correlation effects are taken into account, 
Lee et al. demonstrate a significant suppression in 
the strength of the peak and suggest this change is 
significant for its superconductivity (but to compli- 
cate things, several experimental groups 3 ^— have 
reported the existence, of varying degrees, of anti- 
ferromagnetic fluctuations). Application of DMFT 
to MgFeGe could be very instructive in trying to un- 
derstand differences underlying the different super- 
conducting behavior of these two very similar com- 
pounds. 



III. APPLYING THE MODIFIED 
BECKE-JOHNSON POTENTIAL 

Whereas the observed FSs in several iron- 
pnictides agree rather well with first principles re- 
sults, the agreement has been worse for LiFeAs. Not 
only effective masses, but also sizes of the FSs are 
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FIG. 3: (Color online) Valence radial charge density (in 
a.u.) of Fe in LiFeAs and MgFeGe (top), illustrating 
the extreme similarity not only in the region of the 3d 
peak at 0.7 a.u., but also extending out to 1.4 a.u. The 
tiny difference is quantified by plotting the percentage 
difference (bottom). Both are plotted versus the radial 
distance from the Fe atom center. 



renormalized substantially from calculated results. 
LiFeAs is however a good non-magnetic Fermi liq- 
uid and the effects of interactions as obtained within 
DMFT treatments do not show much more than 
these changes (such as strong satellites, for exam- 
ple). Thus a mean field treatment of the electronic 
structure can be fairly accurate, but it is not the 
one given by the GGA functional. For this reason 
we have investigated the changes given by a recent 
extension of density functional exchange-correlation 
functionals. 

We have applied the mBJ exchange potential, in 
the form known as the TB-mBJ potential (after 
Tran and Blaha's 2 ^ modification), to both LiFeAs 
and MgFeGe. The mBJ potential is an extension of 
the Bccke-Johnson potential^ which itself was for- 
mulated to reproduce the shape of the atomic opti- 
mized effective potential) 41 ' 42 an exact exchange po- 
tential. This goal is approached in a phenomenolog- 
ical manner through the use of orbital kinetic energy 
densities expressed in local-density form. The mBJ 
potential produces rather accurate gaps for a wide 
range of semiconductors and insulators, but as a lo- 
cal density functional it is much more cost-effective 
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FIG. 4: (Color online) Band structures of LiFeAs and 
MgFeGe near the Fermi energy. The differences at T, M, 
and X are discussed in the text. 
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TABLE II: Fe 3d orbital occupation matrix elements 
n mm' °f LiFeAs and MgFeGe, from the LAPW sphere of 
radius R = 2.47 a.u. 



than hybrid or GW calculations^ in which the gap 
corrections are comparable in accuracy. 

The mBJ potential has not however been tested 
on many metals or scmimetals. Singh 4 ^. studied 
LaFeAsO using the mBJ potential and concluded 
that while the 10% smaller bandwidth agreed bet- 
ter with experiment than without mBJ, the slight 
changes in the FSs were within the margins of un- 
certainty when one compares an experimental to a 
relaxed As position. One can note that LaFeAsO 
and the two compounds we study here are semimet- 
als (i.e., "insulators" with slightly overlapping va- 
lence and conduction bands 2 ^) so a functional that 
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FIG. 5: (Color online) Total DOS comparisons between 
the non-mBJ (regular GGA) and mBJ calculations of 
LiFeAs (above) and MgFeGe (below). 



makes favorable corrections for insulators might lead 
to similar improvements for semimctals. 



A. Bands and DOS 

The DOS spectra of LiFeAs and MgFeGe from cal- 
culations with and without the use of mBJ are com- 
pared in Fig. O where it can be seen that the posi- 
tion and intensity of peaks is altered in many energy 
regions. In LiFeAs, the distribution of Fe weight 
is influenced strongly by the mBJ potential in the 
range of —2.5 eV to 1 cV, due to the restructuring 
of some bands, shown in Fig. [6] Applying the mBJ 
potential, the bands that make the electron pock- 
ets surrounding the zone corners move up to create 
a van Hove singularity at ep. Angle-resolved pho- 
toemission spectra (ARPES) of LiFeAs reveal a van 
Hove singularity at the Fermi energy^ this is how- 
ever not born out of the zone-corner electron FSs, 
as is the case in our mBJ calculation, but rather the 
two smaller hole FSs centered at T (see Fig. [7]) . Ac- 
cording to ARPES measurements, these particular 



FSs are quite a bit smaller than those derived from 
a DFT calculation. The fiat band at 0.25 eV along 
the F-Z line maintains its flatness but is shifted to 
slightly below the Fermi energy. This is the same 
band that that crosses the FS to create the largest 
hole pocket at T, so its shifting below ep results 
in the hole FS disappearing. The two (originally) 
smaller hole enclosures increase in size while the elec- 
tron pockets overall become smaller under the mBJ 
potential (see also Fig. [7] for the pre- and post- mBJ 
FSs of LiFeAs). 

The same trends can be observed in MgFeGe. 
Changes due to the mBJ potential occur between 
—2 eV and 1.2 eV, including the formation of a van 
Hove singularity at the Fermi level, again caused by 
the increase in states at Ef due to the shrinking elec- 
tron FSs, shown in Fig. [S] The largest hole pocket 
is significantly reduced in size and the smaller two 
are enlarged. The electron pockets take on a pecu- 
liar shape in which the "waist" is squeezed such that 
the M point is no longer enclosed by the FSs. 



B. Orbital selectivity of the differences 

There is intense interest in the multi-orbital na- 
ture of the iron-pnictides and specifically how much 
each one contributes near the Fermi level. There 
is also a need to understand in more detail orbital- 
dependent effects due to the TB-mBJ exchange po- 
tential. Figs. M and HQ] dissect the DOSs of LiFeAs 
and MgFeGe respectively, into Fe d orbital pro- 
jected densities of states (PDOS). Our GGA re- 
sults are consistent with those of Nekrasov and 
collaborators.— The tig orbitals dominate the Fermi 
level for any exchange-correlation functional for both 
systems. Using mBJ, (i) the d z 2 and d xy states shift 
downward by ~ 0.5 eV, whereas the d x 2_ y 2 move 
upward by a similar amount; (ii) the d xz / yz occu- 
pied bandwidth narrows and the peak near e p moves 
to the Fermi level. Note that it is these d xz / yz or- 
bitals that are responsible for the sharp peak struc- 
ture near or at ep. As for the differences between 
materials, the shift of the d xy states is larger for 
LiFeAs, whose substantial d^-orbital PDOS at ep 
drops as a result by an order of magnitude, from 
0.75 eV" 1 to 0.07 cV" 1 . The d xzjyz van Hove sin- 
gularity makes the value of N(£p) very sensitive to 
details (in exchange-correlation potential, probably 
to structure, etc.). 

Borisenko et aZ.— noted that the smaller T- 
centered hole FSs of LiFeAs derived from DFT are 
not realistic, their being much smaller than observed 
in ARPES data. Two separate studies, by Yin et 
als^ and Ferber et al.,— have shown that using 
DFT+DMFT improves on the DFT band structure, 
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FIG. 7: (Color online) FSs in the Brillouin zone, from 
a) a non-mBJ and b) an mBJ calculation of LiFeAs. T 
is located at the unit center. 



MgFeGe 




FIG. 6: (Color online) Band structures of LiFeAs and 
MgFeGe without (solid, black) and with (dashed, green) 
mBJ implementation. 



yielding FSs more consistent with experiment: the 
outer pocket at T is larger and the inner pockets 
are much smaller. Looking at Fig. it may at first 
glance seem as if the mBJ potential affects the three 
hole pockets in the same way DMFT affects them. 
But, as mentioned before, upon observation of the 
bands crossing the Fermi energy along T-X and Z- 
R, it is evident that quite the opposite happens: the 
largest pocket disappears below the Fermi energy 
and the two smaller pockets become larger. The 
mBJ potential acts upon the MgFeGe system in a 
similar fashion. It decreases the sizes of the outer 
hole cylinder and the two electron pockets, and in- 
creases the size of the smallest hole pocket. Over- 
all, the mBJ-applied calculation generates less two- 
dimensional FSs, and does not improve agreement 
with experiment. 



a) GGA 



b) mBJ 



mm 



FIG. 8: (Color online) FSs in the Brillouin zone, from a) 
a non-mBJ and b) an mBJ calculation of MgFeGe. F is 
located at the unit center. 



C. Changes in orbital occupancies 

In both cases of LiFeAs and MgFeGe, when the 
mBJ potential is added, the charge surrounding the 
Fe atom becomes more anisotropic, concentrating 
more in the direction of the As/Ge atom. Despite 
the visible change in the density contours (not shown 
here) upon application of mBJ, the total 3c? oc- 
cupation within the MT sphere of LiFeAs remains 
unchanged (6.65 for our sphere radius), and that 
of MgFeGe decreases almost imperceptibly (6.64 to 
6.63). Thus there is a only a redistribution between 
3d orbitals. Indeed, as shown in Table ITTA the d z 2 oc- 
cupation increases whereas that of d xz / yz decreases, 
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FIG. 9: (Color online) Fe partial DOSs of LiFeAs with 
and without TB-mBJ. 



FIG. 10: (Color online) Fe PDOSs of MgFeGe without 
mBJ, with TB-mBJ, and with S-mBJ. 



by 7-10% for both compounds. 

While the mBJ formalism is not intrinsically spe- 
cific to gapped materials, the parametrizing of the 
coefficients is done empirically by using insulators 
and semiconductors to match the bandgap. Tran 
and Blaha's justification for using orbitally inde- 
pendent potentials to obtain orbitally dependent 
bandgaps is that gapped systems typically have a 
small overlap between the occupied and unoccupied 
orbitals. This is not the case for metals, however, 
and the mBJ implementation may give bogus results 
if the unoccupied orbital potentials are characteris- 
tically different from the occupied. 



IV. SUMMARY 

In this study we have made a comparison of sev- 
eral aspects of the electronic structures of LiFeAs 
and isoelectronic MgFeGe. A repositioning of a DOS 
peak a few tenths of an electron-volt above the Fermi 



level is the most obvious difference, likely arising due 
to differences in Fe hybridization with Ge compared 
to As. As noted in the original report^ the band 
structures near the Fermi level are very similar. We 
have further established that the Fe total 3d orbital 
occupations for the Fe atom are identical for the two 
compounds, and individual orbital occupations dif- 
fer by very little — interesting but too little to be 
significant. The As-Fe-As and Ge-Fe-Ge angles, al- 
though far from being tetrahedrally regular as is the 
case in the higher T c FeSCs, do not differ signifi- 
cantly in these two compounds. The band structures 
are in fact more similar than in other FeSCs that 
show larger variation but comparable superconduct- 
ing properties. MgFeGe's higher density of states 
at the Fermi energy is significant, but argues that 
MgFeGe should be the better superconductor, rather 
than being nonsuperconducting. 

We have tested the recently extended TB-mBJ 
exchange potential for these compounds, since it 
is an exchange potential that seems to work well 
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LiFeAs 





-2 


-1 





+ 1 


+2 


m \. 








-2 


0.683 


0.000 


0.000 


0.000 


-0.109 


-1 


0.000 


0.565 


0.000 


0.000 


0.000 





0.000 


0.000 


0.763 


0.000 


0.000 


+ 1 


0.000 


0.000 


0.000 


0.554 


0.000 


+2 


-0.109 


0.000 


0.000 


0.000 


0.665 


MgFeGe 


\m' 




-1 





+ 1 


+2 


m \. 


-2 


-2 


0.677 


0.000 


0.000 


0.000 


-0.106 


-1 


0.000 


0.590 


0.000 


0.000 


0.000 





0.000 


0.000 


0.734 


0.000 


0.000 


+ 1 


0.000 


0.000 


0.000 


0.576 


0.000 


+2 


-0.106 


0.000 


0.000 


0.000 


0.658 



TABLE III: Fe 3d density occupation matrix elements 
n mm i of LiFeAs and MgFeGe under the TB-mBJ poten- 
tial. 

for non-magnetic insulators. The differences have 
been quantified, but do not seem to resolve the con- 
flict with experimental (ARPES) observations. Thus 
exchange effects between occupied and unoccupied 
non-magnetic states seem to be treated well with 
TB-mBJ, whereas magnetic fluctuations, which are 
also exchange in origin and are treated rather well 
within DFT+DMFT, are not treated well by the 
TB-mBJ approach. 

Our results have not revealed any clearly essen- 
tial features that leave MgFeGe nonsupcrconduct- 
ing while LiFeAs is an impressive, if not yet well 
understood, superconductor with a relatively high 
T c . We suggest that further exploration in parallel 
of these two compounds, by experiment and theo- 
retically, will provide one of the most promising ap- 
proaches to identifying the superconducting mech- 
anism in LiFeAs and, possibly, in the full class of 
FcSCs. 
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VI. APPENDIX: COMMENTS ON THE 
TB-MBJ POTENTIAL 



Recently, improvements on the original mBJ po- 
tential (TB-mBJ) have been made to refit and 
create new sets of parameters based on the type 
of solid (semiconductor, ionic insulator, noble-gas 
solid, etc.) and also by incorporating more gapped 
systems 4£ As LiFeAs and MgFeGe are, of these 
types, closest in character to semiconductors, we 
adopted the semiconductor parameters and applied 
this new potential, referred to as the S-mBJ po- 
tential, to the DFT GGA calculations on our com- 
pounds. The difference between the LiFeAs calcu- 
lations using TB-mBJ and S-mBJ is negligible and 
hence not shown here, but it is more noticeable in 
MgFeGe. Upon observation of the PDOSs of Fe in 
MgFeGe in Fig. [TUJ it is apparent that the S-mBJ 
system is an intermediary of the non-mBJ and TB- 
mBJ systems. This makes sense since the parame- 
ters of the TB-mBJ potential were determined using 
larger-gapped materials as well. 



Implementing the mBJ potential on 
LiFcAs/MgFcGe increases the occupations of 
Fe and As/Ge but decreases the Li/Mg occupation. 
While the TB-mBJ and S-mBJ methods result 
in very similar band structures for LiFeAs, in the 
case of MgFeGe, more noticeable shifts in energy 
of the electronic states can be observed. The 
parameter c in Eq. 7 of Ref. [46| is determined 
by the fitted parameters A and B (for S-mBJ, 
A = 0.267 and B — 0.656), as well as the averaged 
gradient of the density, g, which differs from system 
to system. For LiFeAs g = 1.39 a.u. -1 , and for 
MgFeGe g = 1.27 a.u. -1 . As g deviates further 
from ~ 1.4 a.u. -1 , so does c sem i from ctb- Iu the 
case of LiFeAs, c som i ~ ctb, but this is not true of 
MgFeGe. 



Y. Fuseya, T. Kariyado, and M. Ogata, J. Phys. Soc. 
Jpn. 78, 023703 (2009). 

2 H. Zhai, F. Wang, and D.-H. Lee, Phys. Rev. B 80, 
064517 (2009). 

3 T. A. Maier, S. Graser, P. J. Hirschfeld, and D. J. 



Scalapino, Phys. Rev. B 83, 100515 (2011). 

4 C. Fang, Y.-L. Wu, R. Thomale, B. A. Bernevig, and 
J. Hu, Phys. Rev. X 1, 011009 (2011). 

5 P. M. R. Brydon, M. Daghofer, C. Timm, and 
J. van den Brink, Phys. Rev. B 83, 060501 (2011). 



9 



6 T. Saito, S. Onari, and H. Kontani, Phys. Rev. B 82, 
144510 (2010). 

7 S. Zhou, G. Kotliar, and Z. Wang, Phys. Rev. B 84, 
140505 (2011). 

8 M. Berciu, I. Elfimov, and G. A. Sawatzky, Phys. Rev. 
B 79, 214507 (2009). 

9 H. Takahashi, Y. Imai, S. Komiya, I. Tsukada, and 

A. Maeda, Phys. Rev. B 84, 132503 (2011). 

10 L. Boeri, O. V. Dolgov, and A. A. Golubov, Phys. 
Rev. Lett. 101, 026403 (2008). 

11 T. Egami, B. Fine, D. Parshall, A. Subedi, and 
D. Singh, Adv. Condens. Matt. Phys. 2010, 164916 

(2010) . 

12 Y. Kamihara, T. Watanabi, M. Hirano, and 
H. Hosono, J. Am. Chem. Soc. 130, 3296 (2008). 

13 S. V. Borisenko, V. B. Zabolotnyy, D. V. Evtushinsky, 
T. K. Kim, I. V. Morozov, A. N. Yaresko, A. A. Ko- 
rdyuk, G. Behr, A. Vasiliev, R. Follath, et al., Phys. 
Rev. Lett. 105, 067002 (2010). 

14 K. Kusakabe and A. Nakanishi, J. Phys. Soc. Jpn. 78, 
124712 (2009). 

15 S. Deng, J. Kohler, and A. Simon, Phys. Rev. B 80, 
214508 (2009). 

16 C. He, Y. Zhang, B. P. Xie, X. F. Wang, L. X. Yang, 

B. Zhou, F. Chen, M. Arita, K. Shimada, H. Na- 
matame, et al., Phys. Rev. Lett. 105, 117002 (2010). 

17 L. Ma, G. F. Chen, D.-X. Yao, J. Zhang, S. Zhang, T.- 
L. Xia, and W. Yu, Phys. Rev. B 83, 132501 (2011). 

18 K. Kitagawa, Y. Mezaki, K. Matsubayashi, Y. Uwa- 
toko, and M. Takigawa, J. Phys. Soc. Jpn. 80, 033705 

(2011) . 

19 D. R. Parker, M. J. Pitcher, P. J. Baker, I. Franke, 
T. Lancaster, S. J. Blundell, and S. J. Clarke, Chem. 
Commun. pp. 2189-2191 (2009). 

20 C. Chu, F. Chen, M. Gooch, A. Guloy, B. Lorenz, 
B. Lv, K. Sasmal, Z. Tang, J. Tapp, and Y. Xue, 
Physica C 469, 326 (2009). 

21 I. Morozov, A. Boltalin, O. Volkova, A. Vasiliev, 
O. Kataeva, U. Stockert, M. Abdel-Hafiez, D. Bom- 
bor, A. Bachmann, L. Harnagea, et al., Cryst. Growth 
Des. 10, 4428 (2010). 

22 X. Liu, S. Matsuishi, S. Fujitsu, and H. Hosono, Phys. 
Rev. B 85, 104403 (2012). 

23 C.-H. Lee, A. Iyo, H. Eisaki, H. Kito, M. T. 
Fernandez-Diaz, T. Ito, K. Kihou, H. Matsuhata, 
M. Braden, and K. Yamada, J. Phys. Soc. Jpn. 77, 
083704 (2008). 

24 D. Johrendt, H. Hosono, R.-D. Hoffmann, and 
R. Pottgen, Z. Kristallogr. 226, 435 (2011). 

25 F. Tran and P. Blaha, Phys. Rev. Lett. 102, 226401 



(2009) . 

26 S. Lebegue, Z. P. Yin, and W. E. Pickett, New J. 
Phys. 11, 025004 (2009). 

27 P. Blaha, K. Schwarz, P. Sorantin, and S. Trickey, 
Comp. Phys. Comm. 59, 399 (1990). 

28 K. Schwarz and P. Blaha, Comput. Mater. Sci. 28, 
259 (2003). 

29 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

30 J. H. Tapp, Z. Tang, B. Lv, K. Sasmal, B. Lorenz, 
P. C. W. Chu, and A. M. Guloy, Phys. Rev. B 78, 
060505 (2008). 

31 D. J. Singh, Phys. Rev. B 78, 094511 (2008). 

32 I. Nekrasov, Z. Pchelkina, and M. Sadovskii, 88, 543 
(2008). 

33 M. Wang, M. Wang, H. Miao, S. V. Carr, D. L. Aber- 
nathy, M. B. Stone, X. C. Wang, L. Xing, C. Q. Jin, 
X. Zhang, et al., ArXiv e-prints (2012), 1208.0909. 

34 G. A. Sawatzky, I. S. Elfimov, J. van den Brink, and 
J. Zaanen, Europhys. Lett. 86, 17006 (2009). 

35 G. Lee, H. S. Ji, Y. Kim, C. Kim, K. Haule, 
G. Kotliar, B. Lee, S. Khim, K. H. Kim, K. S. Kim, 
et al., ArXiv e-prints (2012), 1205.6526. 

36 P. Jeglic, A. Potocnik, M. Klanjsek, M. Bobnar, 
M. Jagodic, K. Koch, H. Rosner, S. Margadonna, 
B. Lv, A. M. Guloy, et al., Phys. Rev. B 81, 140511 

(2010) . 

37 Z. Li, Y. Ooe, X.-C. Wang, Q.-Q. Liu, C.-Q. Jin, 
M. Ichioka, and G. qing Zheng, J. Phys. Soc. Jpn. 
79, 083702 (2010). 

38 L. Ma, J. Zhang, G. F. Chen, and W. Yu, Phys. Rev. 
B 82, 180501 (2010). 

39 A. E. Taylor, M. J. Pitcher, R. A. Ewings, T. G. Per- 
ring, S. J. Clarke, and A. T. Boothroyd, Phys. Rev. 
B 83, 220514 (2011). 

40 A. D. Becke and E. R. Johnson, J. Chem. Phys. 124, 
221101 (2006). 

41 R. T. Sharp and G. K. Horton, Phys. Rev. 90, 317 
(1953). 

42 J. D. Talman and W. F. Shadwick, Phys. Rev. A 14, 
36 (1976). 

43 D. J. Singh, Phys. Rev. B 82, 205102 (2010). 

44 Z. Yin, K. Haule, and G. Kotliar, Nature Mater. 10, 
932 (2011). 

45 J. Ferber, K. Foyevtsova, R. Valentf, and H. O. 
Jeschke, Phys. Rev. B 85, 094505 (2012). 

46 D. Koller, F. Tran, and P. Blaha, Phys. Rev. B 85, 
155109 (2012). 



10 



